Effective screening and the plasmaron bands in Graphene 
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Electron-plasmon coupling in graphene has recently been shown to give rise to a "plasmaron" quasiparticle 
excitation. The strength of this coupling has been predicted to depend on the effective screening, which in turn 
is expected to depend on the dielectric environment of the graphene sheet. Here we compare the strength of 
enviromental screening for graphene on four different substrates by evaluating the separation of the plasmaron 
bands from the hole bands using Angle Resolved PhotoEmission Spectroscopy. Comparison with GqW-KVA 
predictions are used to determine the effective dielectric constant of the underlying substrate layer. We also show 
that plasmaron and electronic properties of graphene can be independently manipulated, an important aspect of 
a possible use in "plasmaronic" devices. 



I. INTRODUCTION 

Bohm and Pines^ first discussed the role of interactions be- 
tween charged quasiparticles and the collective charge den- 
sity oscillations of the electron gas (plasmons) in the renor- 
malization of the electronic band structure. Lundqvist^ later 
predicted that the electrons, holes and plasmons could cou- 
ple to form a new quasi-particle known as a plasmaron. 
These plasmarons have been observed experimentally by 
optical^, tunnelling^ and angle resolved photoemission^ 
spectroscopies. In studying the valence electronic structure of 
graphene, side bands due to plasmaron excitation are accessi- 
ble to experimental investigation in great detail^. The depen- 
dence of these features on electron density confirms quantita- 
tively, by comparison to GqW-RPA theory, that these are in 
fact due to electron-plasmon coupling. The resulting charge 
density waves are expected to oscillate at much higher fre- 
quencies than hole charge density waves^ providing a moti- 
vation for the development of "plasmaronic" devices, which 
merge photonics and electronics. A fundamental requirement 
of these devices is the ability to manipulate the electronic and 
plasmaronic properties independently. Graphene, on account 
of its simple and well understood band structure, provides the 
ideal system to investigate these issues^^iiS. 

The dependence of the plasmaron coupling strength on 
the effective screening of graphene by the underlying sub- 
strate has been investigated theoretically^, however no ex- 
perimental work has been published. Here we present An- 



gle Resolved Photoemission Spectroscopy (ARPES) measure- 
ments of graphene on SiC with four different interface layers: 
(6a/3 X 6V3)i?30° C-SiC (henceforth 6\/3 C-SiC), Au-SiC, 
E-SiC and H-SiC. These are expected, from previous experi- 
mental results, to exhibit considerable differences in the ef- 
fective screening of the graphene. We characterize the ef- 
fect of differing substrates on the plasmaronic bands and show 
that in all cases the spectral function of n-type graphene de- 
pends only on a single, substrate dependent, coupling constant 
aQ regardless of the doping level. Analysis of the measure- 
ments indicate a range of substrate effective dielectric con- 
stants (es ~ 87 to ~ 7.8) indicating that, in addition to many 
other unique properties, graphene is ideally suited to investi- 
gate the role of substrate screening on the band structure. 

In the first hexagonal carbon layer formed by thermal de- 
composition of the SiC(0001)^^~^^ surface, every third atom 
is covalently bonded to the substrate. This creates a (6\/3 x 
6\/3)i?30° C layer that does not show the linear band struc- 
ture characteristic of graphene. The subsequent layers formed 
by further thermal decomposition are decoupled from the sub- 
strate/ underlying 6>/3 C layer and exhibit the characteris- 
tic bandstructure of graphene ^^. Previously hydrogen^"— 
, fluorine ^^ and gold^^ have been shown to intercalate under 
layers of 6>/3 C on SiC, decoupling it from the substrate and 
transforming it into graphene. While the single graphene layer 
on 6>/3 C- SiC is n-type with the Dirac crossing, £^d , at ~ - 
0.45 eV^^ graphene formed by intercalation is p-type with £^d 
at 0.15 eV (Au-SiC), 0.18 eV (H-SiC) and 0.79 eV (E-SiC). 



To investigate the plasmaronic bandstmcture, these interca- 
lated systems were converted to n-type by chemical doping 
with K atoms. 



II. BACKGROUND 

The, n-type, Dirac cone picture of the electron dispersion 
near the Fermi energy, ^p , in graphene is shown in Fig. [T] 
A for the non-interacting, single-particle description. In this 
case the bands are linear and form two cones which meet at a 
single point at the Dirac energy, £^d • This picture does not 
take into account the coupling between the elementary charges 
and the plasmons which leads to the creation of a second plas- 
maronic dispersion^i^a^^*^^. The new picture is schematically 
shown in Fig. [T]B where the hole dispersion, shown in red, is 
intersected by a second (plasmaronic) dispersion, shown in 
black. The Dirac crossing is split into three distinct crossings: 
the hole band crossing at Eq , the plasmaron band crossing 
at E2 and the ring like crossing of the hole and plasmaron 
bands at Ei .The Dirac energy, £^d , is the crossing energy 
of the hole bands in the hypothetical non- interacting picture 
which can not be determined from the ARPES data directly. 
To overcome this in the following we will use the crossing 
of the electron- hole bands in the interacting case, £^0 , as an 
approximation. 

The linear band structure of graphene leads to a linear scal- 
ing of the Random Phase Approximation (RPA) dielectric 
function, and hence the momentum and energy separation of 
the hole and plasmaron bands, with kp and ^ ^6,7,10,21^22 ^ jj^j^ 
separation can be determined directly from ARPES data by a 
line shape analysis of the spectral function in terms of a hole 
band (top) and plasmaron (bottom) band. When taking an en- 
ergy vs. momentum cut (E{kx), at ky = 0) the bands form a 
"diamond" as seen schematically in Fig. [T]C from which the 
scaled energy separation, SE = | £^2 — £^0 1 / 1 ^0 — ^f | , and the 
scaled momentum separation, Sk = |/c+ — /c_|/|/cf|, can be 
determined. 

In the absence of phonons, the plasmaronic band is pre- 
dicted to merge with the hole band, and loses spectral weight 
close to the Fermi level, £^f • This is because the electrons 
only couple strongly to plasmons with the same velocity, a 
condition that cannot be met close to £^f • The plasmon and 
the hole band travel with different velocity, greatly reducing 
the composite particle's lifetime. With phonons, the situation 
is more complicated, since the bare plasmon breaks up into 
two modes^^"^^ . It is expected that this breakup will affect 
the spectral function near E^ but to our knowledge these ef- 
fects have not been calculated to date. 

In the one-particle Green function. Random Phase Approx- 
imation (GqW -RPA ) theory developed by Polini et al.^ and 
Hwang et al32J^ the effective environmental screening of 
the graphene layer, e, is related to a graphene effective cou- 
pling constant, ac, via the relatione = e^/aQhvY ~ 2.2/ aQ. 
For clarity note that e is different from the RPA dynamical 
screening function, e(g,cj), and that vp ~ 1 x 10^ m/s is 
the Fermi velocity. Increasing aQ decreases the RPA dielec- 
tric function and therefore increases the separation of the hole 




FIG. 1. The linear band dispersion in; A, the non-interacting sin- 
gle particle (bare band) description, and B, the case of interactions 
between elementary charges and plasmons. In C a linear cut of 
the Dirac region in B {ky = 0) is shown, hole bands are shown 
in red while plasmaron bands are shown in black. The Dirac en- 
ergy crossing, E-d , in A is split into three crossings by the inter- 
actions with the plasmons; Eq , the hole band crossing, E2 , the 
plasmaron band crossing, and Ei the ring like crossing of the hole 
and plasmaron bands. The scaled momentum and energy widths of 
the diamond region in C are given by 6k = |/c+ — /c-|/|/cf| and 
SE = \E2 — Eo\/\Eo — Ef\ respectively. In D a plot of the graphene 
Brillouin zone, showing the principle points and principle directions 
referenced in the paper, is presented. 



(upper) and plasmaron (lower) bands^. The effective envi- 
ronmental screening of the graphene layer, e, is related to the 
dielectric constant from the substrate, es, and the dielectric 
constant from the vacuum above, evac ^ 1, via the relation 
e = (es + evac)/2 => es ^ 2e — 1. In reality this theory ne- 
glects remote band contributions to screening in the graphene 
layer and the polarizability of the K atoms used to electron 
dope the graphene, and therefore what is found is an upper 
limit to the dielectric constant. 



III. EXPERIMENTAL 

Photoemission spectra were obtained from graphene on 
6\/3 C-SiC(OOOl) samples prepared by the high pressure Ar- 
gon method, of Emtsev et al.— and Ostler et aL^^, and 6\/3 
C-SiC(OOOl) samples^ that were annealed in the presence of 
H, F or Au. Hydrogen^"—, fluorine^ and gold^ have been 
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FIG. 2. Experimental spectral function of graphene on (6 a/3 x 6V^)R30° C-SiC, i, graphene on Au-SiC, ii, graphene on H-SiC, iii, and 
Graphene on F-SiC, iv, in th& K — K direction, A, and the Fermi surface, B. The graphene on (6a/3 x 6A/3)i?30° C-SiC sample is n-doped 
by the substrate (ne ^ 1.2 x 10^^ cm~^) while the graphene on Au-SiC (rih ~ 7 x 10^^ cm~^), graphene on H-SiC (rih ~ 6 x 10^^ cm~^) 
and Graphene on F-SiC (rih ~ 4.5 x 10^^ cm~^) are all p-doped by the substrate. 



shown to intercalate under the 6>/3 C layer, decoupling it from 
the SiC substrate such that the n electrons are free to form the 
linear bands characteristic of graphene. The hydrogen inter- 
calation process follows the method of Riedl et aL^^ with 6>/3 
C-SiC(OOOl) samples annealed at -570 °C for - 75 minutes 
in H2. Treatment of the 6^3 C-SiC(OOOl) with fluorine fol- 
lows the previously published method of Walter et al.^^, where 
the sample is heated to — 200 °C in the presence of a XeF2 
crystal and a sacrificial Mo plate. The sample is then soni- 
cated for several hours in ethanol to remove a molybdenum 
oxide layer formed during the process leaving a fluorine inter- 
calated graphene layer. To intercalate the Au approximately 
10 monolayers of Au atoms was deposited on the 6>/3 C- 
SiC(OOOl) samples which were then annealed to — 600 °C 
for rsj 3 minutes, similar to the approach of Gierz et. al^^^ 

ARPES was performed on all four samples at the Elec- 
tronic Structure Factory endstation (SES-R4000 analyzer) at 
beamline 7 of the Advanced Light Source, Lawrence Berke- 
ley National Laboratory. A photon energy of 95 eV was used 
for photoexcitation with an overall energy resolution of — 25 
meV and momentum resolution of ^ 0.01 A. The graphene 
on 6^3 C-SiC(OOOl), H-SiC and F-SiC intercalated samples 
were transported through air to the experimental chamber and 
annealed to - 500 °C, - 400 °C and - 500 °C respectively 
to remove surface contaminants. The graphene on Au-SiC 
samples were prepared in a connected vacuum chamber prior 
to measurement and annealed to ~ 300 °C. During the mea- 



surements the samples were cooled to ~ 20 K using a liquid 
He-cooled cryostat and the pressure was < 2 xl0~^^ Torr. 

Electron doping was induced in all four samples by de- 
positing small amounts of potassium onto the surface from an 
SAES getter source. Charge is transferred from the randomly 
located potassium atoms to the carbon atoms in the graphene 
lattice, leading to an increase in the electron doping of the 
graphene. Addition of progressively more potassium atoms 
allows for any doping level between the starting doping level 
and a limit of~ 6 x 10^^ e/cm^. Above this limit the density 
of potassium atoms cluster into a disordered version of the K 
2x2 structure^^ which provides little doping to the graphene 
and is sufficient to significantly scatter the outgoing electrons 
leading to a broadening of the graphene n bands. 



IV. RESULTS AND DISCUSSION 

A. Spectral function of clean samples 

The experimental spectral functions for the four samples (i - 
iv) considered are shown in Fig. [2]A, for the i^ — i^ direction. 
Renormalization of the bands by phonons in each case is indi- 
cated by the increase in spectral intensity, and a corresponding 
"kink" in the bands between Ep and ~ 200 meV^^. For the 
n doped spectra (graphene on 6>/3 C-SiC sample. Fig. [2] A.i) 
it was previously assumed that the bands below and above the 
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FIG. 3. Experimental spectral function of graphene on Au-SiC, i, graphene on (6a/3 x 6V^)R30° C-SiC, ii, graphene on F-SiC, iii, and 
graphene on H-SiC, iv, in the K — K direction, A, and the K — T direction, B. As opposed to[T]the order is now in terms of increaseing ac 
. All samples are n-doped to ne ~ 6 x 10^^ cm~^ by the addition of K atoms to the surface. All of the plots are presented with the energy 
scaled to the Dirac energy (^scaled = Eb/Eq) and the momentum scaled to the Fermi vector (/cscaied = A;/|A;f|). The results of the energy 
line profile peak fitting to the electron hole bands and the plasmaron bands are shown in red and black respectively, while the A:/|/cf | = line 
profile and fitted peaks are shown in C. The separation of the hole and plasmaron bands increases from left to right and in all cases two bands 
(hole and plasmaron) are observed with the the hole(plasmaron) band showing greater intensity above (below) the Dirac crossings. 



Dirac crossings are offset by an energy shift, SE^^-^^. This 
offset causes a "kink" at the Dirac crossing ( ~ -0.5 eV). A 
significant variation in the position of the Dirac energy below 
(Fig. [2] A.i) and above (Fig. [2]A.ii, A.iii and A.iv) the Fermi 
energy is evident for p and n doping of the graphene respec- 
tively, with hole doping increasing from left to right. 

The doping variation is also clear from the size of the ex- 
perimental Fermi surfaces presented in Fig. [2]B.i - B.iv. The 
Dirac cone picture described above (Fig. [T]) implies that the 
Fermi surface in each case should be circular, this is only the 
case for the Au-SiC (Fig. [2]B.ii) and H-SiC (Fig. [2]B.iii) sam- 
ples. However, the Dirac picture of the bandstructure is only 
valid close to the K point, i.e. for the small Fermi surfaces 
in Fig. [2]B.ii and B.iii. As the momentum values increase 



away from the K point the three-fold symmetry of the Bril- 
louin zone around the K point results in the trigonally warped 
Fermi surfaces observed in Fig. [2]B.i and B.iv. 

The anisotropy in the Fermi surface intensity has also been 
studied^ ^'^^ and is related to the chiral nature of the electron 
states in graphene on account of the two sublattices. The ex- 
tinguishing of intensity in the negative (positive for n-doped 
samples) ky direction is strongly affected by symmetry break- 
ing, as is seen in bi-layer graphene samples^^'^^. Therefore the 
lack of intensity in the negative (positive for n-doped samples) 
ky direction in the current samples rules out strong symmetry 
breaking due to the substrate^ ^'^^ and the associated forma- 
tion of a gap from causing the observed offset of the upper 
and lower bands as was proposed previously^. This approxi- 
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FIG. 4. Constant energy ( Ei ), i, and constant momentum (k = 
0), ii, intensity profiles as a function of K doping for Graphene on 
Au-SiC, A, Graphene on (6^ x 6V3)i?30° C-SiC, B, Graphene 
on F-SiC, C, and Graphene on H-SiC, D, in the K — T direction. 
All of the plots are presented with the energy scaled to the Dirac 
energy (^scaled = Eb/Eq) and the momentum scaled to the Fermi 
vector (fecaied = /c/l/crD-The momentum ( 6k ) and Energy ( 6E 
) separations in i and ii are constant except close to the zero doping 
level, where the vanishing density of states act to converge the hole 
and plasmaron bands. 



mation has been shown to be valid at the high photon energies 
employed here^~— . 



B. Influence of electron doping on the spectral function 

Let us recall that the separation between the hole and plas- 
maron bands is constant as a function of chemical doping 
when normalised to £^d and /cp . Here we analyze this sep- 
aration for four different substrates as a function of doping. 
In Fig. Owe present the experimental spectral functions from 
each of the samples, doped to ne ~ 6 x 10^^ cm~^ by the ad- 
dition of potassium to the surface in the iiT — i^ (Fig. [3] A) 
and K — r directions( Fig. [3]B), with energy and momen- 
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FIG. 5. Peak amplitudes and Gaussian widths from the peaks 
shown in Fig. [3] for Graphene on Au-SiC, A, Graphene on (6 a/3 x 
6y3)i?30° C-SiC, B, Graphene on F-SiC, C, and Graphene on H- 
SiC, D, in the i^ — F direction. In all cases we observe a similar 
increase (reduction) in the amplitude (Gaussian line width) of the 
hole band peaks, shown in red. The reverse is true for the plasmaron 
band peaks, shown in black. The fitted positions from Fig. [5] clearly 
show that for ky / kp < (ky / kp > 0) the positions are below 



(above) the Dirac crossing. 



tum scaled to ^d and k^ respectively. Line shape analysis 
is performed by non-linear least squares fitting of the spec- 
tral function at fixed wave vector, k, for which an example is 
given in Fig|3]C (k=0) for each of the samples. The results of 
peak fitting to the hole (upper) and plasmaron (lower) bands 
are represented by the red and black dotted lines respectively. 
The constant energy (binding energy ( Eb ) = £^d ) and con- 
stant momentum (k = 0) intensity maps from spectra similar 
to FigUB are plotted as a function of potassium doping in Fig. 
|4]i and ii. The normalized width of the spectra, corresponding 
to the energy and momentum width of the "diamond" region, 
in Fig. [4] are remarkably constant except at low doping where 
the scaling relationship is expected to break down and phonon 
interactions are expected. As this scaling relation is unique to 
Coulombic electron-electron interaction effects this provides 
strong support for the model. 

Representative Energy Distribution Curves (/c/|/cf| = 0) 
obtained from the spectral functions in Fig[3]B are presented 
in FigOC. Lorentzian- Gaussian peaks are fitted to the pro- 
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FIG. 6. The energy separations, SE , (in black) and momentum 
separations, Sk , (in red) of the "diamond" region as a function of 
graphene effective dielectric constant, ao , as determined from one 
particle Green's function calculations within the Random Phase Ap- 
proximation (GqW -RPA y . The values for ac that correspond to 
the experimental widths shown in Table U are indicated by dashed 
lines. 



files indicating that two bands are observed in all cases. The 
Lorentzian linewidth for each peak is fixed at '^ 0.1, while all 
other parameters are fitted. The separation of the two peaks 
is obvious in the fluorine and hydrogen intercalated samples, 
Fig. [3]C iii and iv, however it is clear that a second component 
is needed to arrive at a reasonable description in the case of 
the Au intercalated and 6>/3 C samples. The positions of the 
peaks determined from these fits are overlay ed on Fig 15] A and 
B in red (hole band) and black (plasmaron band) and provide 
a good description in all four cases. At high binding energy 
only the plasmaron band is observed, the hole band emerges 
from the plasmaron band close to the lower Dirac crossing, 
and the two bands are separated. After passing through all 
three Dirac crossings the plasmaron band curves back up to 
meet the hole band, with the separation of the bands at the 
Dirac crossing, Ei, being different for each of the samples. 
The amplitude (circles) and Gaussian linewidth (crosses) of 
the fitted peaks are shown in Fig. \5\ We observe a decrease 
in the amplitude, and corresponding increase in the linewidth, 
of the plasmaron band with increasing momentum. The re- 
verse is true for the hole band. This behaviour is predicted by 
the GqW -RPA plasmaron theory of Polini et al. and Hwang 
et aliiiSiSiiS^ and provides further support to the existence of 
plasmaron bands in graphene on all four substrates. 

There is strong debate about the source of the deviation of 
the upper and lower bands in graphene on 6\/3 C-SiC, with 
both the plasmaron bands-^ and the existence of a symmetry 
breaking-induced gap^i^Si^ proposed. Our work shows that 
the variation of the shape of the bands in Fig. 3 can be ac- 
counted for purely by electron-electron interactions (In any 
case, symmetry breaking would be expected to be very weak 
for graphene on Au, H, or F). Qi et al3 have presented a 



defect model for the structure of the 6>/3 C layer based on de- 
fects observed in STM data, and another defect based model 
was presented by Kim et al^. In the Qi model^S a signifi- 
cant deviation from the structure of defect free graphene is 
proposed and therefore a large potential barrier would need 
to be overcome to transform the 6^/3 C layer into defect free 
graphene. This is in contrast to the relative ease with which 
atomic layers are intercalated under the 6>/3 C layer (eg. Au^^« 
, F^^ and H^^ presented here) resulting in essentially defect 
free graphene. 

Qi et s& and Kim et al^^ propose a structure for the 6\/3 C 
layer, based on STM measurements, and use this to calculate 
a band structure for both the 6\/3 C layer and the graphene 
on the 6\/3 C layer. Contributions to the electronic structure 
from the 6>/3 C layer close to the Dirac energy are used to ex- 
plain the apparent gap in the ARPES spectra of graphene on 
6>/3 C-SiC. The experimental spectral function of the buffer 
layer has been investigated previously ^^ and none of the de- 
tailed ARPES spectral features predicted by the models of 
Qi et al^ or Kim et al^ were observed. Importantly, neither 
model correctly predicts the insulating behaviour of the 6>/3 
C layer, with no states above -0.4 eV binding energy observed 
experimentally^^. Defect models cannot, therefore, explain 
the experimental observations of either the 6^/3 C layer or 
graphene on 6\/3 C- SiC, and the deviation from linear be- 
haviour at £^D is due to the formation of plasmaronic bands in 
all four of the samples presented here, including graphene on 
6^3 C-SiC. 



C. Effective screening analysis 

In Fig. [6] the variation of the width and height of the dia- 
mond region ( SE and Sk respectively) with graphene effec- 
tive coupling constant, ac , as determined from GqW -RPA 
theory is shown. These values were determined by generating 
a theoretical spectral functions for a range of ac values and 
fitting peaks to determine the height and width of the diamond 
region. The experimental widths, heights and the correspond- 
ing coupling constant (derived from these by comparison with 
theory) are indicated by dashed lines and recorded in Table [D 
As described in section |II1 the effective coupling constant can 
be used to determine the dielectric constant, e, and an upper 
limit to the substrate effective dielectric constant, es, which 
are also listed in Table H The metallic substrate, Au, provides 
the highest effective screening (es ^87), an order of magni- 
tude larger than F (es ^10) or H (es ^7.8). Interestingly the 
supposedly electronically decoupled 6>/3 C layer (es ^43) 
provides a relatively high effective screening. This can be at- 
tributed to the presence of dangling Si bonds, which are quite 
polarizabla^. 



D. Comparison to GqW -RPA theory 

In order to provide a detailed comparison between theory 
and experiment we present in Fig. [7] the calculated spectral 
functions on the basis of the GqW -RPA theory of Polini et 
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FIG. 7. One particle Green's function calculations within the Random Phase Approximation ( GoW -RPA ^ of the spectral function of 
graphene in the K — K (A) and K — F (B) directions for the effective coupling constants (ac = e^ /ehvp ^ 2.2/e ) found experimentally. 
Overlayed are the experimental fitted positions from FigOfor Graphene on Au-SiC, i, Graphene on (6\/3 x 6v^)^30° C-SiC, ii, Graphene on 
F-SiC, iii, and Graphene on H-SiC, iv. Above the lower Dirac crossing, E2 , good agreement between the experiment and theory is observed, 
however below E2 the hole (red) and plasmaron (black) bands quickly hybridize, deviating from the experimental prediction. The intensity at 
high binding energy is also over-estimated by the theory (see Fig. [3])- 



TABLE I. Effective screening constants for Graphene determined 
from experimental spectral functions: Energy ( 6E ) and momentum 
( Sk ) separation of the hole (upper) and plasmaron (lower) bands, 
graphene effective coupling constant ( ac — e^ /ehvF ^ 2.2/e ), 
the graphene effective dielectric constant (e) and the substrate effec- 
tive dielectric constant (es ~ 2e — 1) 



Substrate 


6E 


6k 


ac e 


es 


Au-SiC 


0.12±0.04 


0.09±0.02 


0.05±0.01 44±9 


87±18 


eVsc-sic 


0.21±0.02 


0.16±0.01 


0.1±0.04 22±8 


43±16 


F-SiC 


0.40±0.05 


0.29=b0.03 


0.4±0.05 5.5±0.7 


10±1.3 


H-SiC 


0.49±0.02 


0.34±0.01 


0.5±0.03 4.4±0.3 


7.8±0.5 



al.- in the K — K ( A ) and i^ — F ( B ) directions us- 
ing the experimentally determined alpha values. Overlayed 
are the experimental fitted hole (red) and plasmaron (black) 
bands obtained for Graphene on Au-SiC, i, Graphene on Qy/S 
C-SiC, ii, Graphene on F-SiC,iii, and Graphene on H-SiC, iv. 
In all cases good agreement is found between theory and ex- 
periment above the high energy Dirac crossing, E2 . Below 
E2 the experimental hole and plasmaron bands quickly merge 
in contrast to the prediction of the GqW -RPA theory. The 
experimental intensity (see Fig. [3]) in this high energy region 
is also poorly described by the theory indicating that further 



theoretical work needs to be undertaken to accurately describe 
the spectral function below E2 . 

The data presented indicate that the separation of the hole 
and plasmaronic bands can be modified by the use of var- 
ious substrates, which alter the effective dielectric coupling 
in the graphene. Doping of graphene is possible without ex- 
erting an influence on the band separation by the addition of 
potassium to the surface, thereby allowing the plasmaronic 
and electronic band structures to be separately manipulated. 



V. CONCLUSION 

The shape of the spectral function of graphene on 
SiC(OOOl) and the presence of side bands to the linear bands 
near the K point of the Brillouin zone has recently been in- 
terpreted in terms of the formation of coupled hole-plasmon 
quasiparticles (plasmarons). Here we study graphene, on four 
different interface structures (the (6\/3 x 6^/3)R30° recon- 
struction, and the SiC - graphene interface intercalated with 
gold, hydrogen or fluorine), using angle-resolved photoemis- 
sion. Our data show that similar side bands are found in all 
cases as revealed by energy and momentum distribution curve 



line shape analysis. These show the scaling of the separation 
with the energy of the Dirac point and the Fermi wave vector. 
The separation of the plasmaronic side bands and the bands 
unaffected by hole-plasmon coupling is found to strongly vary 
with interface structure. We use GqW -RPA calculations to 
determine an upper limit to the effective dielectric constant of 
the underlying substrate interfacial layer; this is found to vary 
from es ~7.8 for hydrogen to es ~87 for Au intercalation. 
Graphene is thus an ideal candidate for investigating the ef- 
fective screening in the context of hole-plasmon interactions. 
We also show that plasmaronic and electronic properties of 
graphene can be separately manipulated. 
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